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Abstract 

We study the discrete and gauge symmetries of Quantum Electrodynamics at 

finite temperature within the the real-time formalism. The gauge invariance 

of the complete generating functional leads to the finite temperature Ward 

identities. These Ward identities relate the eight vertex functions to the 

elements of the self-energy matrix. Combining the relations obtained from 

the Z2 and the gauge symmetries of the theory we find that only one out 

of eight longitudinal vertex functions is independent. As a consequence of 

the Ward identities it is shown that some elements of the vertex function are 

singular when the photon momentum goes to zero. 
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I. INTRODUCTION 



Quantum field theory at finite temperature (QFFT) has been a subject of great current 
interest for a variety of physical contexts 0. Two main formulations of QFFT can be 
distinguished. In the imaginary-time formalism (ITF) |2|] the continuous energy variables 
k are replaced by the discrete Matsubara frequencies u n , and loop integrals over k are 
transformed into discrete sums. To obtain a graph with real external energies requires a 
non-trivial analytic continuation. On the other hand, the real-time formalism (RTF) allows 
the computation of dynamical quantities directly as functions of continuous real frequencies. 
In addition, it is possible to formulate the theory in a covariant way [[j]. The price to 
paid is that the Feynman rules are now more complicated that those at zero temperature. 
The number of degrees of freedom is doubled and hence the propagator acquires a matrix 
structure. The RTF has been discussed using a canonical operator formulation || and also 
using the path integral method ||. Here, we shall mainly consider the diagrammatic method 
based in the path integral representation of the generating functional. 

One would expect physical quantities to be independent of the formalism used to calcu- 
late them. This, however has been an issue that has generated a longstanding controversy. 
For example, the amputated n-point functions at finite temperature appear to lead to dif- 
ferent results when calculated in the two formalism's. Nevertheless, it has been shown that 
within the real-time formalism one can consider a sum of graphs, motivated by causality ar- 
guments, which at least for the two- and three-point functions agree with the corresponding 
analytically continued imaginary-time results Thus, whereas single graphs may differ 

in the two formalism's, there exist quantities that agree when calculated by either method. 

In this work we are interested in studying the discrete and gauge symmetries of QED 
at finite temperature within the real-time formalism. It is commonly believed that gauge 
invariance is not affected by the temperature ||. There exist, of course, some previous 
studies on the subject ||||; in particular, in the case of gauge theories at high temperature 
it has been shown that the hard thermal loops and the tree amplitudes obey the same Ward 
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identities [Hj. That Ward identities may be carried over to finite temperature without any 
modification, except for the fact that the matrix elements depend on the temperature, is 
probably obvious in the ITF. However in the RTF, where there are eight different vertex, the 
structure of the Ward identities is far from obvious. To our knowledge an explicit derivation 
of the Ward identities in the RTF has not appeared in the literature. The purpose of this 
article is to provide a systematic derivation based on the path integral method. In particular, 
we consider the Ward identities that relate the 2- and 3- point Green functions. 

Section |T| is devoted to the study of the effective two component thermal formulation 
of QED within the path integral formalism's. Along the paper we stress the fact that the 
thermal information in the action is contained in the convergence factor e that is inserted 
through boundary conditions. We also discuss the discrete symmetries of the effective QED 
action. In section |T| the subject of the gauge invariance of the theory is examined. Ward 
identities are obtained by demanding the complete generating functional be gauge invariant. 
When applied to the photon propagator, these identities imply that any component of the 
photon polarization tensor is transverse to any order of perturbation theory. Furthermore, 
the Ward identities impose several relations among the vertex functions and the elements 
of the self-energy matrix. We exploit these relations to prove that only one out of eight 
components of the longitudinal part of the vertex functions is independent. We also show 
that the Ward identities at finite temperature demand that some of the vertex functions 
have to be singular when the photon momentum goes to zero. Finally, the general results 
are illustrated by a one-loop calculation. 

II. THERMAL QED 

Let us consider a gas of electrons and photons at finite temperature. The system is 
described by the QED Lagrangian 

C = --F^F^ + + ieA^ - - |(<9 • A) 2 , (2.1) 
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which includes the gauge fixing term. The path integral representation of the generating 
functional can be written as 

Z [J„, 77, 77] = J [DA^LhpLhp] exp {ij^drdx (c + fa + # + J V A V ) } , (2.2) 

where the integration over r is along a monotonically descendent contour C in the complex 
time-plane. The contour has end points T; and T; — z/3, with /? denoting the inverse of the 
temperature. The fields in Q2.2|) satisfy the boundary conditions A u {ji) = A v {ji — if}) , and 
if(Ti) = —e^if{Ti — if}), where /i is the chemical potential. Except for the gauge-fixing and 
the source terms the action in (|2.2j ) is gauge invariant. Indeed, the first three terms in L are 
invariant under the gauge transformation 

Ay -> A„ + d„A , i) -> exp {-ieA} $ , (2.3) 

with the time derivative taken along C. 

Within the RTF the curve C is chosen in such a way that includes the real time axis. A 
particular contour that has been extensively used is one of the form C = C\ © C 2 © C 3 © C 4 . 
The segment C\ starts at Tj = —T and runs along the real time axis up to +T. C<i continues 
along the imaginary axis from +T to +T — iaf}, C3 goes from +T — iaf} to — T — iaf}, and 
C4 runs from — T — iaf} to — T — At some stage, the limit T — ► 00 is taken to produce 
the standard RTF. The contour parameter a can take any value within the interval (0, 1). 
Here we shall adopt the symmetrical value a = 1/2. Of course, physical results should not 
depend on a. 

Following the method of Niemi and Semenoff 0] Z can be recast as an integral over fields 
evaluated at real time, but with the number of degrees of freedom effectively doubled. The 
fields lying on C\ (type-1 field) are referred as the physical field, while the (type-2) fields 
associated with the section C2 of the contour are named as (thermal) ghost fields. In terms 
of these fields the generating functional takes the form jTIJ : 



Z [J£, Va , fj a ] = J [DA£DiJ> a D$ a ] exp {iS eff } , (2.4) 
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where the effective action 



S eff = J (^t 1 ! 2 ] + i'aVa + Va^Pa + JvcA^) dx , 



(2.5) 



has been written in terms of the finite-temperature Lagrangian £[1,2] = + Ca + £int, 
with 



5 (0) 



a6 



^6, 



^6 > 



(2.6) 



(0) 



In these equations a summation over the repeated index a = 1, 2 is understood, while 5^ 



and 



£>(0) 



denote the matrix elements of the bare thermal propagators. At the tree level, 
there are only two kinds of vertices, the type-1 and the type-2 vertices, that differ by a sign: 
r^ n = — r^ 22 = 7 M . The physical Green functions can be expressed as a sum of Feynman 
diagrams, with the type-1 fields appearing on the external legs; the type-2 field only occur 
in the internal lines of the diagrams. 

In momentum space the inverse of the fermion propagator can be worked out as 



S(°)(p)l 1 = y>-m)a 3 + ie\V- 1 (p ) 



(2.7) 



where 03 is the diagonal Pauli matrix and 

cos (p p — e (po ) e^/ 2 sin ip p ' 



V = 



e(po)e ^ 2 sm(f p 



cos (f p = 



6(p )e x / 4 + e(-p )e- x / 4 



cos ip p 



^ e x/2 + e -x/2 



(2.8) 



with x — (3{p — [/,). For the photon the result is 



-1 



-k 2 gp, + (1 - OM 5 * 



0-3 ~ icguv 



U-^fco) 



(2.9) 



with 



U = 



cosh^i. sinh^i 



sinh 9 k cosh # fc 



cosh 2 6 k = 



1 - exp-^l fc ol 



(2.10) 



We have kept e finite in order to define the theory properly At zero temperature e 
is included as a convergence factor, so the path integral in Minkowski space is well defined. 
At finite temperature the procedure becomes essential, not only for convergence, but also 
to keep trace of the temperature dependence, which as seen from the previous expression 
appears only in the e-terms. Because of these contributions, the action is non-local in the 
time variable. In fact, the substitution of ( p.7|) in (|2.6|) leads to the following expression for 
the fermion part of the action in configuration space 

£</, = J dx dyodx^ a (x ,x) [5(x - y ) - m) a 3 + ieM(x - y )] ab ip b (yo, x) , (2.11) 



where 



dk 

M ab (x -y )= I ^ [V-\h 



Similarly, for the gauge-field we obtain 



La = J dxodyodx A^(xq, x) x 
1 



-ik (x -y ) 



ab 



[5(x - yo ) {QtxvU + (f - l)d^d u }a z + ieg /1 „N(x - y )] ab A u b (y 0: x) 



(2.12) 



(2.13) 



with 



N afe (x -y ) 



dk 
27 



-ik (x -y ) 



(2.14) 



In the zero temperature limit N = M = 15 (xo — yo), where I is the identity matrix. Then, 
the physical and ghost fields uncouple and the time-locality of the action is restored. 

In reference it was pointed out that the scalar theory at finite temperature has a 
Z 2 symmetry under the interchange of type-1 and type-2 fields. Here we exploit the Z 2 
symmetry for thermal QED to reduce the number of independent 2- and 3-point vertex 
functions. As usual we introduce the generating functional of one particle irreducible graphs 



v]> vl/ 4^ 

^ai -"-a 



W [fj a , Va, J'] ~ J (Va^ + « W + J£A$ 



(2.15) 



where W [r) a ,T] a , J£\ = —UriZ [r) a ,r) a , Jff] is the connected generating functional and 

A 4 ' A Aa h* ' 



5 J, a " ' 8 A 



8W _ 8T 
8fj a ' <5^ a ~~ ^ 



f = £ = (2.16) 

In Eqs. ( |2.15| ) and ( |2.16| ), \l/ a , \I/ a , and A% refer to the classical fields, but for simplicity we 



have denoted them with the same notation than the one used for the fields appearing in the 
path integral. It is convenient to write the fields as two component spinors 

A^= , ¥= , * = (Vi ^2 • (2.17) 

In terms of them the Z 2 symmetry for T is realized as follows: 

A**] = -r*[**5-, 5-**, cr 1 A /1 ] , (2.18) 

where a\ is the Pauli matrix and 

/ o ieM 2 \ 

a=\ . (2.19) 

As is well known, the fermion 2-point function 8 2 T / 'S^Sip can be related to the inverse 
of the exact propagator. The explicit relation in momentum space is 

{2irf5{p-p')iS^{p) = f dzdyjV*-^ T (2.20) 



Taking into account that S 1 (p) = S^°\p) — SQo), the Z 2 symmetry (|2.18|) imply that 



the elements of the fermion self-energy matrix S are related as follows: Sn(p) = — £ 22 (p), 
^12 (f) = e ^2i(p)- An additional relation among the elements of of S is obtained from the 
circling relation ||: 

En(p) + e^''/ 2 S 12 (p) + e^"/ 2 £ 21 (p) + £ 22 (p) = . (2.21) 



The previous relations taken together imply that only one component of the matrix self- 
energy is independent. Therefore, for example, knowing Eu is enough to determine the other 
components of S. In particular, for the retarded self-energy, E#(p) = En(p) + e~ /3p0//2 Ei 2 (p), 
we have : 

ReE R (p) = ReSu(p) 

bn£ fl (p) = e(po) sec(2p p ) ImE u (p) = -ie~M 2 csc(2y? p ) E 12 (p) . (2.22) 

In practice it is simpler to calculate ImE# in terms of the E12, which is a purely imaginary 
quantity . For this reason, and also for future comparison with the result for the 3-point 



function (see Eq. ( p.!6|) ), we have written the second equality for the imaginary part of E#. 

We now turn the attention on the proper vertex function. Within the RTF there are 
eight connected three point functions T^ bc (q,p), which are defined by 

ie (W-p- dl^,,,) = J dM V^ Jf _^__ . (2 . 23) 

The first index corresponds to the photon with momentum q, while the second and the third 
indices refer to the incoming fermion with momentum p and the outgoing fermion with 
momentum p + q, respectively. Not all of these complex functions are independent. The Z 2 
symmetry of the theory ( |2.18| ) leads to the following relations: 

1^222 (<?,p) = - \Piu(q,p)]* , 

^ 22 (q,p) = -[TU<l,P)Y ■ (2-24) 

The circling equation of Ref. || adapted to the vertex function of QED introduces a further 
relation among the various T^ bc , namely 

J- 111 -r e 211 ' 121 ' 1 221 

+ e ^/ 2 r^ 12 + e^r^ + e^/ 2 r^ 22 + r& 2 = o , (2.25) 



where r = po + qo- Equations (|2.24j) and (|2.25|) represent five complex constrictions and 



consequently, one is left with only three independent 3-point functions. It is convenient to 
choose them in the following way |J: 

rkfop) = r£u + e-^/ 2 r^ 21 + e^ 2 T^ l2 + e^ 2 r? 22 , 

T^( q ,p) = 1% + e~^ 2 n u + e^' 2 T» l2 + e^T^ , 

r» m (q, P ) = r? u + e-**/*^ + e-^r^ + e-^ 2 r% 21 . (2.26) 

These are the retarded functions with, respectively, tx, £2, and t 3 as the largest time, their 
complex conjugate being the corresponding advanced functions. In the next section we shall 
see how the Ward identities lead to a further reduction in the number of independent 3-point 
functions. 



III. GAUGE INVARIANCE AND WARD IDENTITIES 

As stated in section O, at finite temperature the Lagrangian (without gauge-fixing) with 
support in the complex contour C is gauge invariant. However, the original theory has 
been replaced by an effective one with two components, and we must ask about the gauge- 
invariance properties of this effective theory, as well as of their consequences. To investigate 
this point, let us consider the following independent gauge transformations on the type-1 
and type-2 fields: 

i[) a -> exp(-ieA a )^a , (3.1) 

where A a (x) (a = 1, 2) are arbitrary functions of x, that in what follows we assume to 
be infinitesimal. Under such transformation, the finite temperature effective action S e ff 
in ( |2.5| ) is not invariant; the origin is threefold: The same as in the zero temperature 
case, the gauge fixing term and the source term are changed; in addition, the temperature 
dependent contributions, contained in the e-terms, are also gauge dependent. Accordingly, 
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it is convenient to decompose the variation of the effective action into two parts and write 
5S e ff = 5 S + 5 St- 5S arises from the variation of the gauge fixing and source terms, and 
using (|2.6|) , ( |2.11| ) and ( |2.13| ) it can be worked out as 



ss 



^ ay ^ ay ^ a 



J dx -£e a (d • A a ) □ + J^dft - ie (j] a ip a - tpaVc 



A„ 



(3.2) 



with ei = 1 , e 2 = — 1. The quantity 5 St is obtained from the e-terms in ( [2.1 1| ) and ( (2.13|) . 
The result is 



5S 



T 



* ay * ay -"-a 



■>< j r/.;W.r (l (/// n <j - A a (x )-^ [N ab (x - y )A b M (y )] 



+ie [A a (x Q ) - A b (y )} ip a (x , x)M ab (x - y )i/; b (yo, x) 



(3.3) 



Although S e ff varies with the gauge transformations, the physical consequences of the 
theory, expressed in terms of Green's functions, should not depend on the gauge. Thus, the 
generating functional Z must be gauge invariant and this non-trivial requirement leads to 
the Ward identities. Combining the previous results with Eq. ( |2.4| ) the gauge invariance of 
Z implies 



5 S 



5 5 5 
Srja' 5fja 5 Jo 1 



5S 7 



5 5 5 
Srja 5fja 5 JL 1 



Z[J^,VyV] = 0. 



(3.4) 



By means of Eqs. (|2.15| ), (|2.16|) and using the results in (|3.2|) and ( |3.3|) we convert ( [3.4|) into 
a condition for T: 



SAZix) 



$i>a(x) 



Ha(x) 



+iee J dyo 



dtp (y ) 

itp a (xo)M ab (x - yo)ip b (y ) + M ab (x - y ) ~ b , °, - U' 



dVa(yo) 



b, xq y ) 



d 

- ie j d yog^I [ N ab(x - y )A^(y , x)] = 



(3.5) 



where we used that the gauge functions A a (x) are arbitrary, and no summation over the 
index a is implied. 

Equation ( |3.5|) expresses the general content of the Ward identities in the real-time for- 
malism. Repeated differentiation of it, at A% = ip a = if> a = 0, generates relations between the 
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one particle irreducible Green functions that are consequence of the gauge invariance of the 
theory. First we differentiate ( |3.5| ) with respect to A^{x) and using the fact that ST /SA^SA U 
is the inverse of the full photon propagator D^, the following result (in momentum space) 
is obtained 



g M D MI/ 1 (q) = -^a 3 q u q 2 - ieq u U ^go 



(3.6) 



This relation for the full inverse photon propagator holds to all orders in perturbation theory. 
The effect of radiative corrections is to add the negative of the self-energy to the inverse 
of the bare propagator [D°J _1 . From the expression in (2^9) it is easy to check that (515) is 



satisfied by [D° ] 1 . Consequently we get 



qm^q) = 0. 



(3.7) 



In a recent work, the transversality of the polarization tensor at finite temperature has been 

corroborated by an explicitly one-loop calculation Jl2]| . As we proved here, this property 

remains valid to all orders in perturbation theory. 

As a second application of the identity (|3.5|), we act on it with 8 2 /5ip(x)5^(y) and put 
= ip a = ip a = 0. Using the definitions of the proper vertex function ( |2.23| ) and the full 

inverse fermion propagator fl2.20| ), we arrive to 



>ab 



Sj{p + q) - it V x (p + go 



5 n 



Sab(p)-ie \V-\ P0 ) 



ab 



. (3.8) 



These are the Ward identities that, in the RTF, relate the components of the full electron 
propagator to the vertex functions. Equation ( |3.8j ) represents eight relations, one for each 
element of T^ bc . They include explicit temperature contributions due to the presence of the 
e-terms. At zero temperature, the matrix V(po) reduces to the identity, the e dependent 
terms cancel, and ( ^8|) decompose in two independent relations, one for the type-1 fields 
and another for the type-2 fields. 

If we substitute the inverse bare propagator ( ^7|) , Eq. (|3~8l ) simplifies to 



% T abM,P) = q^CaSabSc 



(3.9) 
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which are trivially satisfied, because the only vertices that do not vanish at the tree level 
are T^q = — r 222 = 7 M - According to this result, it is convenient to write 



KbMp) = T^A&^a + A^ bc (q,p) 



(3.10) 



where A^ fec (g,p) represents the radiative corrections to the vertex function. In terms of 
the fermion self-energy matrix the inverse of the complete fermion propagator can 



be written as S 1 (p) = S^°'(j>) 
implies that 



S(p). This relation combined with Eqs. (|3~8l ) - (|3.10|) 



Explicitly, these relations are 



(3.11) 



q^ A lll(q,p) 


= -[£ n (p + g) 


-S u (p)] 


% A ii2(q,p) 


= -E 2l (p + q), 




M-i2i(g>p) 


= £12 (p) , 




9^122(9^) 


= 0, 




Q^ A 2ii(Q,P) 


= 0, 




Q^212(l,P) 


= £ 21 (p) , 




Qti A 22l(l,P) 


= -Ei 2 (p + q), 




g M A£ 22 (g,p) 


= - [S 22 (p + g) 


- S 22 (p)] 



(3.12) 

Several important results can be read from the above identities. First, the vertex func- 
tions with equal thermal fermion indices but a different photon index (i.e. Aj 1 22 and A 211 ) are 
transverse. Next, in the limit — > the diagonal vertex functions A^ 11 (0,p) and A 222 (0,p) 
may be determined as follows in terms of the diagonal components of the self-energy: 

0£ii(p) 



A? n (0,p) 
A£ 22 (0,p) 



dZ 22 (p) 
dpu 



(3.13) 
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Furthermore, comparing the second (third) line with the sixth (seven) line in ( ft. 12] ) we see 
that the relations 

A&i(<Z,p) = -Kai{q,P + q), (3-14) 

are satisfied by the longitudinal components of the corresponding elements of A^ bc . Notice 
that the fermion momentum in the right side is shifted from p to p + q. The four vertex 
functions in the last equations are imaginary since S 12 and S 21 in Eqs fl3.12| ) are imaginary 
quantities. 



Equation (|3.12| ) implies that, either the non-diagonal elements of S vanish or that the 



vertex functions with different values for the fermion thermal index have a singularity at 
= 0. Since we know that the self-energy is in general non-diagonal, we conclude that the 

four vertex function in (|3.14|) are singular in the g M — > limit. Later on, this fact will be 

illustrated within a one loop-calculation. 

The foregoing relations can be applied to simplify the longitudinal part of the effective 

casual vertex defined in ( |2.26|) . First, if we combine the circling equation ( |2.25|) with ( 2.24|) 



and ( |3.14j ) we can derive the following relation 



ImA1 n (q,p) = ie~ M2 [csc(2^ p+(? )A? 21 (g,p + q) - csc(2^ p )A? 21 (g,p)] , (3.15) 

for the longitudinal components. Using these results we then find that the retarded vertex 
functions in Eq. (|2.26|) can be written as 



ReA^p) = ReA^p) = ReA^p) = ReAt n (q,p) , 

ImA^(q,p) = X -e-M 2 [csc(2^ p+ ,)A5 t 21 (g,p + q) + csc(2^)A^ 21 (g,p)] , 

ImA^ R2 (q,p) = -ImA^(q,p) = -e^ /2 [csc(2^ p+(? )A^ 21 (g,p + q) - csc(2^ p )A^ 21 (g,p)] . 

(3.16) 

These equations are very useful because they show that we need only know one (complex) 
of eight vertex functions A a b c to determine the retarded vertex function of QED. In fact, 
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according to ( |3.16|) we can choose to calculate Re Ai U (q,p) and Im Aj t 21 (g,p) to determine 



the other components of A^ bc as well as the three retarded vertex functions. Although, 
Ai2i(q,p) is singular for g M — > 0, we notice that this singularity is cancelled for A^ R2 (q,p) 
and Aft 3 (q,p). Consequently, only the imaginary part of A R1 (q,p) presents a singularity at 

Referring back to Eqs. (|2.22|) we see that the real and imaginary parts of E# are obtained 
from i?eSn and Si 2 respectively. Here, we have derived an analogous result for the retarded 
vertex functions. Indeed, in (|3.16|) the real and imaginary parts of the retarded vertex 
functions are expressed in terms of ReAm and A121 respectively. However, it should be 
stressed that whereas the result for the fermion self-energy holds in any theory, Eq. (|3.16|) 
ensues from the gauge invariance of the theory. 



Finally, combining ( |2.26| ) with ( |3.12|) the Ward identities for the retarded vertex functions 
read 

qv A m(q,p) = ~ £U(p + g) - ^r(p)} , 

% A m(q,p) = q^m(q,p) = - [^r(p + q) - ^r(p)} , (3.17) 

where Ha = Y,* R is the advanced self-energy. Thus, we notice that the vertex function 
with the largest time corresponding to the out-going fermion is only related to the retarded 
self-energy. On the contrary, if the photon vertex has the largest time, then the Ward 
identity relates the vertex to a combination of the advanced and retarded self-energies. 
These equations are consistent with the fact that only A^(g,p) is singular at = 0. 

Now, we explicitly verify that Ward identities are satisfied to the one-loop order. To this 
order, the contribution to the vertex function, represented by diagram (a) in Fig. ??, is 



-^A^ bc (q,p) = -0'e) 3 e a e 6 e, 



r d 4 k 1 


D (oy 


/ (2tt)M 





Z(khM°J(P + <l-kh^sl J( 



(3.18) 



The fermion and photon propagators are the bare propagators in Eqs. ( [2.7D and ( |2.9p 
respectively. 
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Similarly, the first contribution to the self-energy is shown in diagram (b) of Fig. ?? 
and reads 



- iT, ab (q,p) = -(ie) 2 e a e b 



r d 4 k 1 


D (oy 


/ (2tt)M 





pa 

6a 



{khJS^J {{> - k)-f p 



(3.19) 




p+q 



(b) 




p-k 



(a) The vertex A^ bc (q,p) to the one-loop order, (b) Fermion self-energy to the one-loop 
order. 



The bare fermion propagator in (|2.7| ) satisfies the relation 
S$t(p)U-m) = eJ ab + o{e> e( ^- m2) 



(3.20) 



^ ' (p 2 — m 2 ) 2 + e 2 j 

Here we can safely take the limit e — > 0, obtaining well defined results. These are all 



that we need. If we contract the vertex in ( 3.18Q with g M , making the replacement $ = 
d + i> ~ ft — m ) ~ (J> ~ ft — m ) an d using Eqs. ( |3.20| ) and ( |3.19| ) we immediately obtain the 
result given in Eq. ( |3.11| ). Hence, the Ward identities in ( |3.11 and |3.12| ) are indeed verified 
at the one loop-order. 

Finally, we turn the attention to the origin of the singular behavior of the vertex functions 
A-abdl, p)i b 7^ c in the limit — ► 0. Let us consider for example the function A 12 i(q,p), 
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an easy algebraic manipulation of Eq. (|3.18 ) shows that the loop integral contains, among 
others, the factors 

sin 2 (p p+q _ k A*(p + q - k)A(p - k) - cos 2 ip p+q ^ k A(p + q - k)A*(p - k) (3.21) 

where A(p) = l/(p 2 — m 2 + ie) is the Feynman propagator. Then, it is obvious that for 
any value q^ ^ the integral is free of pinch singularities, or A*(fc)A(fc) terms. However, 
in the q^ — > limit the integral develops a pinch singularity. This is the origin, at the 
one-loop level, of the singular behavior of the vertex function. Nevertheless, it must be 
stressed that we have proved that the singular behavior should remain valid to all orders in 
perturbation theory. A detailed calculation of the one-loop vertex functions A abc (q,p) with 
the corresponding analysis of the structure of the singularity at g M = will be presented 
elsewhere |13| . 



IV. CONCLUSIONS 

We have studied the discrete and gauge symmetries of Quantum Electrodynamics at 
finite temperature within the the real-time formalism. In the path integral representation of 
the generating functional the thermal information is inserted through boundary conditions, 
which are taken into account in the e factor. The presence of this convergence factor leads 
to an action that is non-local in time, besides being gauge dependent. By demanding the 
generating functional to be gauge independent we obtain the Ward identities, which relate 
the eight vertex functions to the elements of the self-energy matrix. Combining the relations 
obtained from the Z 2 symmetry and the gauge symmetry we find that (for the longitudinal 
part) only one out of eight vertex functions is independent. In addition, the retarded vertex 
functions obey relations similar to those known for the retarded self-energy. 

As a consequence of the Ward identities, the vertex functions A^ bc with b ^ c are found 
to be singular when the photon momentum goes to zero. The zero-momentum limit of 
Feynman amplitudes at finite temperature has generated much discussion. In particular, 
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it is known that for the thermal self-energy calculated in perturbation theory, the limits 
Po — > and \p\ — *■ are not interchangeable. It has been argued that the nonanalyticity of 
the self-energy at p^ = should be expected on physical grounds. In particular Weldon [TJj 



presented various methods to prove that the self-energy for scalar bosons has a branch cut 
in p and consequently the po — > and \p\ — > limits do not commute. On the other hand, 
various attempts have been made to improve the usual perturbation expansion in such a way 



that the self energy becomes analytic at p^ = [15|]. Although our results refer to the 3-point 
function of QED, they may be relevant in the understanding of the zero-momentum limit of 
Feynman amplitudes at finite temperature. The singularity of the QED vertex function at 
zero momentum occurs because of the Ward identities, and therefore is not a consequence of 
the approximation made in a perturbative calculation. However, it should be stressed that 
in this paper the singularity in the QED vertex functions is predicted based on the Ward 
identities, whereas this kind of singular behavior is known to exist even for theories without 
gauge invariance. 
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